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Abstract 

Mg-based alloys have recently attracted major interest in view of their potential 
applications in the aerospace, aircraft and automotive industries. Here, we show 
that the effects of Y and Zn atoms on their compressibility can be reliably esti- 
mated by a simple thermodynamical model deduced by means of first-principles 
calculations based on density functional theory that just appeared. 
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1. Introduction 



Due to the high importance in modern technology, the effects of sohd solution 

strengthening in Mg alloys has been recently investigated (e.g., sec Refs [1-3]) 
using various kinds of solute atoms. After the observation [4] that the hardness 
of Mg-Zn alloys increases upon increasing the Zn content, additional investi- 
gations [5-7] revealed the strengthening effect of Zn in Mg-Rarc Earths (RE) 
solid solutions. The latter solutions have attracted major attention in view of 
their high specific strength at both room and elevated temperatures and of their 
excellent creep resistance [1, 3, 8]. In particular, Y is considered to be one of 
the most important RE elements to improve mechanical properties of Mg-alloys 
at high temperatures [9]. 

Very recently [10] first-principles plane-wave pseudo-potential calculations based 
on density functional theory (DFT) appeared which studied the effects of Y and 
Zn atoms on elastic properties of Mg solid solutions. The bulk modulus B, shear 
modulus, Young's modulus and Poisson ratio were also derived [10] by Voigt- 
Reuss approximation, indicating a strong dependence of the elastic modulus on 
concentration of the solute atoms. Here, focusing solely on the bulk modulus 
(or its inverse, i.e., the compressibility k, k = 1/B), we shall show that values 
comparable to those obtained by the aforementioned first-principles calculations 
can be determined on the basis of a thermodynamical model [13, 14] (for a re- 
cent review, see Ref. [15]). Note that such a procedure has already provided 
successful results for the compressibilities of ionic solid solutions, either of NaCl- 
or CsCl-type structure (see Refs [16] and [17] respectively), as well as for the 
transition-metal carbides and nitrides alloys Zra;Nbi_a;C and Zr^Nbi-^N [18]. 

2. The procedure to estimate the compressibility of a solid solution 
in terms of the compressibilities of the end members 

Let us denote Vi and V2 the corresponding molar volumes, i.e. Vi = Nvi and 

V2 = Nv2 (where N stands for Avogadro's number) of two pure end members 
1 and 2 respectively. A "defect volume" [19, 20] v''' is defined as the variation 
of the volume V\, if one "molecule" of type (1) is replaced by one "molecule" 
of type (2). It is then evident that the addition of one "molecule" of type (2) 
to a crystal containing N "molecules" of type (1) will increase its volume by 
v'^ + v\. Assuming that u** is independent of composition, the volume Vn+u of 
a crystal containing A'' "molecules" of type (1) and n "molecules" of type (2) 
can be written as: 

Vn+u = [1 + {n/N)\V^ + nv'^ (1) 

The compressibility n of this solid solution can be obtained by differentiating 
Eq. (1) with respect to pressure which gives [19]: 

kVn+u = KiVi + {n/N) [k'^Nv'^ + KiVi] (2) 

where k'^ denotes the compressibility of the volume f'', defined as 

k'^ = -{1/v'^) X {dv'^/dP)T. (3) 
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To a first approximation, the "defect-volurae" u'' can be estimated from: 

= {V2 ~Vi)/N (4) 

Thus, if Vn+ti can be determined in general from Eq. ([T]) (or approximately 
from Eq. (4)), the compressibility k can be found from Eq. ([2]) in case that a 
reliable estimate of k'^ can be made. 

Along this direction, we employ the thermodynamical model, termed cBVl 
model, for the formation and migration of the defects in solids [13-15]. This 
model has been checked for defect processes in a variety of solids [19] as well as 
in cases of activation of defects in complex ionic materials where upon gradually 
increasing the pressure (P) a critical pressure is reached at which an electric 
signal is emitted well before fracture thus providing an explanation for the de- 
tection of signals before major earthquakes [21, 22]. 

According to the aforementioned thermodynamical model, the defect Gibbs en- 
ergy (7* is interconnected with the bulk properties of the solid through the rela- 
tion = c*i3J7 where B stands, as mentioned, for the isothermal bulk modulus 
(= 1/k), the mean volume per atom and d is dimensionless quantity. (The 
superscript i refers to the defect process under consideration, e.g. defect for- 
mation, defect migration and self-diffusion activation). By differentiating this 
relation in respect to pressureP, we find the defect volume [= {dg^ /dP)T]- 
The compressibility k'*'* defined by k''^*[= —{d£nv^ / dP)T], is given by [19, 23]: 



K 



{l/B) - {d^B/dP')/[{dB/dP)T - 1] (5) 

We now assume that the validity of Eq. holds also for the compressibility 
k"^ involved in Eq. ([2]), i.e., 

k'^^ki- {d^Bi/dP^)/[{dBi/dP)T - 1] (6) 

where the subscript "1" in the quantities at the right side denotes that they 
refer to the pure end member ([1]). Since, in general, the quantities dBi/dP 
and d^Bi/dP'^, can be roughly estimated from the modified Born model, Eq. 
(O can provide an estimate of k"* and therefrom-on the basis of Eq. ^ the 
compressibility k of the solid solution can be found at various concentrations. 
This procedure has been applied in the former cases [16-18]. In the present 
case, however, we do not apply this procedure, because it is not certain that the 
modified Born model can satisfactorily describe the Mg solid solutions under 
consideration. Thus, we apply the following alternative procedure. 
An inspection of Fig. 3(b) of Ref. [10] verifies that the volume V varies almost 
linearly with respect to different concentrations of Y and Zn, thus agreeing with 
the essence of Eq. ([T]). Hence, to a first approximation we can write 

V={1- x)Vi + XV2 (7) 

the differentiation of which in respect to pressure leads to: 

l + .(t-l) 



B = B, 



(8) 
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This relation enables the estimation of the bulk modulus of the solid solution, 
for various concentrations, in terms of the bulk moduli Bi and B2 of the pure 
end members. The application of Eq. ([5]) leads to values that strikingly coincide 
with those computed in Ref. [10] (see their figure 5(a)) on the basis of first- 
principle calculations. 

It should be emphasized that the success of the thermodynamical model [13, 14, 
19] adopted here, has been also checked very recently in several tens of other 
mixed systems (see Ref. [24] along with its improvement in Ref. [25]). 

3. Conclusion 

Using a thermodynamical model which interconnects the point defect parame- 
ters in solids in terms of the bulk properties, we obtained a relation, i.e. Eq. (|8]), 
which allows the estimation of the bulk modulus of solid solutions in terms of 
the bulk moduli of the end members. This gives results that strikingly coincide 
with those deduced very recently for Mg-based alloys by means of first-principle 
calculations based on density functional theory. 
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